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Abstract 

A modified Woods-Saxon potential model is proposed for a unified description of the entrance 
channel fusion barrier and the fission barrier of fusion-fission reactions based on the Skyrme energy- 
density functional approach. The fusion excitation functions of 120 reactions have been systemat- 
ically studied. The fusion (capture) cross sections are well described with the calculated potential 
and an empirical barrier distribution. Incorporating a statistical model (HIVAP code) for describ- 
ing the decay of the compound nucleus, the evaporation residue (and fission) cross sections of 51 
fusion-fission reactions have been systematically investigated. Optimal values of some key param- 
eters of the HIVAP code are obtained based on the experimental data of these reactions. The 
experimental data are reasonably well reproduced by the calculated results. The upper and lower 
confidence limits of the systematic errors of the calculated results are given. 



Electronic address: wangning@gxnu.edu.cn 



I. INTRODUCTION 



15|. Theoretical 



The production of superheavy nuclei as evaporation residues in fusion reactions is a field 
of very intense studies in the recent decades 1,2,3,4,^^ So far^ the superheavy elements 
Z = 110 ~ 116 and 118 have been synthesized [7|, [sl, l9|, [lO|, lul, lli 
support for these very time consuming experiments is vital in choosing the optimum target- 
projectile-energy combinations, and for the estimation of cross sections and identification of 
evaporation residues. A self-consistent microscopic dynamics model is still not yet available 
for practical studies of the whole fusion process from the capture to the decay of the heavy 
compound nuclei. Therefore, in the practical calculation of the evaporation residue cross 
section, the reaction process leading to the synthesis of superheavy nuclei is divided into 
two or three steps. Firstly, the projectile is captured by the target and a dinuclear system 
is formed which then evolves into the compound nucleus, and finally, the compound nucleus 
loses its excitation energy mainly by emission of particles and 7-ray and goes to its ground 
state. The simplified version of the evaporation residue cross section is given by 



(Ecm. ) -PcN (Ec.ra. ) W^sur (-^c 



(1) 



Here, (Tcap, -Pcn and Wsut are the capture cross section for the transition of the colliding 
nuclei over the entrance channel Coulomb barrier, the probability of the compound nucleus 
formation after the capture and the survival probability of the excited compound nucleus, 
respectively. There are several unsolved questions in each component of the right side of 
Eq.(l) which leave a certain margin of uncertainty in the estimates of the evaporation residue 
cross section [l^. In addition, there could be several parameters in the practical models 
which are hardly unambiguously determined by a very limited number of measured evap- 
oration residue cross sections of superheavy nuclei. For example, t^e calculated formation 
probability Pcn of the compound nuclei for reaction ^^Fe+^'^^Pb in [3] is about two orders of 



magnitude larger than that obtained in Ref. |4j], both of the models can, however, reproduce 
the measured evaporation residue cross sections satisfactorily. Therefore, it is necessary to 
test and determine the interaction and parameters adopted in each component of Eq.(l) 
individually. 

To study the three components in Eq.(l) individually, we first investigate the influence of 
the fission and quasi-fission on the fusion-fission reactions. It is generally thought that for 
systems with the compound-nuclear charge number Zcn smaller than about 60, the fission 



barrier is high enough to make fission an improbable decay mode at incident energies close 



to the fusion barrier 



18| . Thus for these reactions, uer — o"fus — o"cap holds at near-barrier 



energies. To see it more clearly, we present a schematic figure (Fig. 1(a)). The horizontal and 
vertical axis denote the compound- nuclear charge number Zcn and the mass asymmetry of 
the reaction system rj = {A2 — A1)/ {A2 + A1), respectively. Here, Ai and A2 are the projectile 
and target masses. The fusion reactions in region I have Pcn — W^ur — 1 as discussed before. 
There are quite a large number of experimental data of evaporation residue cross sections for 
the reactions in region I accumulated in recent decades, which makes it possible to establish 
a reliable model for systematic description of the capture process without the infiuence of 
fission and quasi-fission. For heavier compound systems the fission increases rapidly with the 
Zq^/Aci^} and the angular momentum. For sufficiently asymmetric systems with Zcn well 
below 100 (systems in region II of Fig. 1(a)), and at energies close to the fusion barrier, it is 
generally recognized that (Xfus = crER + crpF- Here the aj^s, c"er and crpF are the cross sections 
for fusion, evaporation residue and fission, respectively. For systems in region II, it is thought 
that the quasi-fission barrier is high enough and thus Pcn — 1- The available experimental 
data of the evaporation residue cross sections for reactions in region II are less than those 
in region I, but they seem to be much enough for a systematic investigation to test and 
determine some key parameters of a statistical model for calculating the survival probability 
WsuT, combining the model for describing the capture cross sections, without the infiuence 
of the quasi-fission. In addition, the measured fusion cross sections for reactions in region 
II can further test the theoretical model for calculating o"cap- For Zqn larger than about 
100 (systems in region HI of Fig. 1(a)), quasi-fission occurs. Thus in the calculation of the 
evaporation residue cross sections for these reactions, the infiuence of quasi-fission should be 
taken into account (Pcn < !)• In Fig. 1(b) we illustrate this point more clearly. We show the 
contour plot of the quasi-fission barrier heights of the reactions with reaction parters along 
the /3-stability line. Here, the height of the quasi-fission barrier is empirically estimated by 
the depth of the pocket of the entrance channel capture potential obtained with a modified 
Woods-Saxon potential which will be discussed in Sect. II. The height of the quasi-fission 
barrier decreases rapidly with the increase of the compound-nuclear charge, especially for 
symmetric target-projectile combination. For reactions with the same 2'cn, those with more 
asymmetric target-projectile combinations have higher quasi-fission barriers. It is expected 
that for sufficiently asymmetric systems the fusion probability Pcn is approximately equal 
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FIG. 1: (Color online) (a) A schematic figure for different types of fusion reactions. The horizontal 
and vertical axis denote the compound-nuclear charge number Zcn and the mass asymmetry of 
the reaction system r/ = {A2 — Ai)/{A2 + ^1), respectively, (b) Contour plot of the quasi-fission 
barrier heights obtained with a modified Woods-Saxon potential (which will be introduced in the 
next section) for reactions with nuclei along the /3-stability line. 

to one as mentioned above. If both cTcap and W^sur can be predicted reliably, this would help 
to understand the dynamics of fusion and quasi-fission. 

Based on above discussions, the emphasis of this paper is put on study of such fusion- 
fission reactions in which the quasi-fission is not important. To study this kind of reactions 
we employ a modified Woods-Saxon potential model based on the Skyrme energy density 



functiona 
posed in 



together with the extended Thomas-Fermi approach. This model was first pro- 
and a large number of fusion reactions have been described satisfactorily with 



the entrance channel potential. The potential between nuclei around the touching point can 
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or calculation of Wgur, the well known standard sta- 
I20I ) is used. Then, the evaporation residue cross 



be accurately evaluated with numerical algorithm . Unfortunately, it is not so convenient 
for any practical application because one needs to evaluate numerically the microscopic den- 
sities of the interacting nuclei, the derivatives of these densities and the integrals. It is better 
to find an analytical expression for the potential. In this work we will present an analytical 
modified Woods-Saxon (MWS) form for the potential based on the numerical results. With 
the analytical MWS potential, both the fusion barrier and the fission barrier of a reaction 
system will be consistently studied, 
tistical model (with HIVAP code 5 
sections of a series of fusion-fission reactions will be investigated for a systematic test of the 
model and refining the parameters. 

II. MODIFIED WOODS-SAXON POTENTIAL AND SOME PARAMETERS 

OF HIVAP CODE 

In this section, we first introduce an empirical nucleus-nucleus potential based on the Skyrme 
energy density functional within the extended Thomas- Fermi approach. Then, the statistical 
model HIVAP is briefiy introduced and the infiuence of some key parameters is studied. 
Finally, a number of calculated results are compared with experimental data. 

A. Modified Woods-ScLxon Potential and Fusion Cross Section 

The nucleus-nucleus interaction potential reads as: 

V{R) = Vm{R) + Vc{R). (2) 

Here, V^v and Vc are the nuclear and Coulomb interactions, respectively. We take Vc{R) = 
e^ZiZ2/ R-, and the nuclear interaction to be of Woods-Saxon form with five parameters 
determined by fitting the entrance channel potentials obtained with the Skyrme energy 
density functional within the extended Thomas- Fermi (up to second order in h 2l|) approach 
proposed in 

Vq 

^^^^^ = l + exp[(/?-/?o)H' 



with 
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^1/3^1/3 

Vo = u,[l + K{h + h)] I ' (4) 

^1 + ^2 



TABLE I: Parameters of the potential. 
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and 



Rn 



l/3x 



(5) 



/i = (A^i — Zi)/Ai and I2 = {N2 — Z2)/A2 in Eq.(4) are the isospin asymmetries of projectile 
and target nuclei, respectively. 

By varying the five free parameters uq, tq, c and a of the modified Woods-Saxon 
(MWS) potential, we minimize the relative deviation between the fusion barrier height 
obtained with the Skyrme energy-density functional with SkM* 23|] force and the barrier 
height of the MWS potential obtained with Eq.(2). The corresponding optimal values of 
these parameters are obtained at the minimum of the relative deviation. In this work, 66996 
reactions with Z1Z2 < 3000 were used to determine the parameters of the modified Woods- 
Saxon potential. The obtained optimal values of the parameters are listed in Table I. Here 
we also list the potential depth parameters Uq and k, proposed in 22(] for comparison. In 22 1 
the nuclear interaction is taken as a Gaussian form and the potential parameters are also 
determined by the Skyrme interaction SkM*. We find that the potential depth parameters 
obtained with the two approaches are close to each other. 

With the modified Woods-Saxon potential together with the proposed empirical fusion 
barrier distribution in [l^, the fusion cross sections and the mean barrier heights of a large 



number of reactions can be reproduced well 



IS 
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261 . For the reader's convenience, the 



empirical barrier distribution is briefiy introduced here. We assume the barrier distribution 
function D{B) to be a superposition of two Gaussian functions Di{B) and D2{B), 



Di{B) 



V7 



exp 



-7- 



(26i)2 



(6) 



and 



D2{B) 



exp 



{B - B2f 

{2h2f 



(7) 
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with 



(8) 



B2 = B, + 62, 



(9) 



B, 



(10) 



b2 = -iBo-B,). 



(11) 



Here Bq is the barrier height from the modified Woods-Saxon potential. The effective barrier 
height is Be = fBo with the reduction factor / = 0.926. The quantity 7 in Di{B) is a factor 
which empirically takes into account the structure effects and has a value larger or equal 
to 0.5. For the fusion reactions with neutron-shell open nuclei but near the /3-stability line 
and for the fusion reactions at energies near and above the barriers we set 7 = 1. For the 
reactions with neutron-shell closed nuclei or neutron-rich nuclei an empirical formula for the 
7 values was proposed in 19|]. For a more convenient discussion, we introduce the inverse of 



7 as an enhancement factor (7 = 1/7. The larger the value of g is, the larger the fusion cross 
section at sub-barrier energies is. From the discussions in 19|], we learn that for the reactions 



with neutron-shell closed nuclei we have < g < 1 while for the reactions with neutron-rich 
nuclei 1 < g < 2. With the proposed empirical barrier distribution, and the fusion radius 
Rfus and the curvature of the barrier hu obtained with the modified Woods-Saxon potential, 
the fusion excitation function (or the capture excitation function of reactions in region III 
of Fig. 1(a)) can be obtained (details in Refs. [19|, l24| ) 



(Ecm.) = min[cri(Ec .m. ) ) C^avr (^c.m.)], 



(12) 



with 



and 



DliB) aJ::\E,,^„B)dB, 



D^{B) + D2{B) 



crJ::'{E..r..,B)dB. 



(13) 



(14) 
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FIG. 2: (Color online) The average deviations X\og ^ total of 120 fusion reactions. Ai and A2 
denote the projectile and target masses, respectively. 



Where, cr^^""^ denotes Wong's formula 27 1 for penetrating an one- dimensional parabolic 



barrier, 

~27l 



Wong/ TP TJ \ ^^}us , A , 

CTf^, n^c.m., Bo) = In (1 + exp 



f. [Ecm. — Bq 

njjj 



(15) 

with the center-of-mass energy -Ec.m.- Bq, Rfus and fuj are the barrier height, radius and 
curvature, respectively. The influence of angular momentum in the entrance channel has 
already been taken into account in Wong's formula with the assumptions that the barrier 
position Rfus and the barrier curvature hu do not change with angular momentum. 

We have calculated the fusion (capture) excitation functions of 120 fusion reactions at 
energies near and above the barrier (with g = 1) and their average deviations X\og from the 
experimental data defined as 

^log = ^Y1 [log(ai,(K)) - log(ae..p(K))]' . (16) 

n=l 

Here m denotes the number of energy-points of experimental data, and ath{En) and aexp{En) 
are the calculated and experimental fusion (capture) cross sections at the center-of-mass 
energy En {En > Bq), respectively. The calculated results for X\og are shown in Fig. 2. The 
average deviations of about 70% systems in Xiog are less than 0.005, with which we can 
estimate the systematic error of this approach for the description of the fusion (capture) 
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FIG. 3: (Color online) The fusion excitation functions of a series of reactions with ^^O bombarding 
on medium mass targets. The incident energies are normalized by the mean barrier heights B^. 
The fusion cross sections of these reactions are shown with linear and logarithmic scale in (a) 
and (b), respectively. The solid circles and crosses denote the experimental data of reactions with 
neutron-shell closed nuclei and with neutron-shell open nuclei, respectively. The solid curve denotes 
the calculation result with g = 1- The error bars in (a) are estimated with 18% of the fusion cross 
sections. The upper and lower limits of the cross sections in (b) are obtained with g = 2 and g — 0, 
respectively. 

cross sections at energies near and above the barriers. A series of fusion reactions with ^^O 
bombarding on medium mass targets such as ^"^^"^^^Sm are studied with this approach, and 
the fusion excitation functions of these reactions are shown in Fig.3. The energy scale has 
been normalized by the mean barrier height -Bm calculated with the proposed method by 
setting g = 1 26|]. The scattered symbols denote the experimental data. The solid curve 
denote the calculated results with g = I. The error bars in Fig.3(a) are estimated by 18%. 
In Fig.3 (b), we notice that nearly all of the experimental data of sub-barrier energies are 
scattered in the region < < 2 as we defined in the proposed approach. The fusion 
cross sections (solid circles) of the reactions with neutron-shell closed nuclei at sub-barrier 
energies are systematically lower than the calculated results with g = 1 which is consistent 
with our discussion mentioned previously. With g — and g = 2 we estimate the lower and 
upper limits of the fusion (capture) cross sections at sub-barrier energies respectively. 



B. Fission Barrier and Level Density Parameter in Evaporation Calculations 



The calculations of the survival probabilities W^ur of the compound nuclei were performed 
with the statistical evaporation code called HIVAP which uses standard evaporation theory 
and takes into account the competition of 7-ray, neutron, proton, a-particle emission with 
fission using angular-momentum and shape- dependent two-Fermi-gas-model level density 
formulaj^. Although it is a standard statistical model for describing the de-excitation 
process, one has to reconsider some parameters adopted for describing a wide range of 
fusion-fission reactions. The sensitive parameters involved are primarily fission barriers and 
level density parameters. 

In the standard HIVAP code, the fission barrier at zero angular momentum is calculated 

by 



B 



f 



oMac 



s. 



(17) 



The macroscopic 



ic barrier B^'^ is usually described with a liquid-drop model refined by Cohen 
and Swiatecki 28|, Sierk 29|], and Dahlinger et al. 30[. The shell correction S is calculated 
from the difference of the experimental mass and the liquid-drop mass, S = M^^p — Mi^d. 
In this code, the liquid-drop mass is calculated with the parameter set proposed by Myers 



and Swiatecki in 1967 



31| , and the Mexp is in fact taken from the mass table of Moller-Nix 



32( 1 which was obtained with the finite range droplet model and has an rms deviation of 



only 0.656 MeV for 2149 measured masses of nuclei 33|. In the present work, we calculate 
the macroscopic fission barriers with the proposed modified Woods-Saxon (MWS) potential 
model in which the parameters of MWS potential are obtained based on the Skyrme energy 
density functional. The value of B^^"^ is empirically estimated by the depth of the potential 
pocket, as shown as an example in Fig.4. This figure is for the 102N0 (formed in reaction 
^'^Ca+^°^Pb) fissioning into two ^IfSb. The obtained barrier is 1.74 MeV. The corresponding 



data from refs. 



281, 129|, 130[ are 1.44, 1.02 and 1.19 MeV, respectively. The barrier for ^^^Pu 



from our method and from refs. [281. 129l. I30l| are 4.16, 5.17, 3.95 and 4.13 MeV, respectively. 
The deviations between our calculated results and the results of liquid-drop models for heavy 
nuclei are in a permitting region. For medium mass nuclei, our results are in agreement with 
those of refs. 3J, l35|| in which the reduction of the liquid-drop barriers was discussed. 

We know that the nuclear shapes during fission are more elongated than during fusion. 
In this empirical approach, the neck and elongation of the system at fission configuration 
can not be described well in the sudden approximation. We concentrate on the height of the 
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FIG. 4: (Color online) The macroscopic fission barrier B^^^ for 102^0 fissioning into two ^liSb 
obtained with the modified Woods-Saxon potential. 

fission barrier in this method. We will systematically investigate 51 fusion-fission reactions 
with the fission barriers obtained with four different models (MWS potential model, Cohen- 



Swiatecki's[28|], Sierk's 
the following paragraph 



291 ] and Dahlinger's SOj methods). The results will be discussed in 



In this code, the level density is [20| 

.2\ 3/2 



E*) = - l^-j (2J + l)a'/'Uf exp[2{aUjy/% 



Uj 



E* -Er{J). 



(18) 



(19) 



Here Er{J) is the yrast energy of either the equilibrium configuration (light-particle and 
7-emission) or the saddle-point configuration (fission) and reads 



Er{J) = J{J+l)h^/2I, 



(20) 
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in which / is the moment of inertia. The level density parameter a is obtained from [18| as 

a = d[l + fiE*)S/E*], (21) 

with 36 1 

fiE*) = l-expi- E*/E,) (22) 

with the shell damping energy Ed being 18.5MeV [2^. In the standard HIVAP code, the 
smooth, shell-independent level-density parameter reads 

a = 0.04543 rlA + 0.1355 rlA^^^Bs + 0.1426 r^A^/^S^, (23) 

which takes into account the volume, surface and curvature dependence of the single-particle 



level density at the Fermi sur: 
defined in the droplet model 



ace. Bs and Bk denote the surface and curvature factors 



3^]. For evaporation channels we set Bs = Bk = 1. For 
the fission channel, the values of Bs and B^ are tabulated as a function of the fissility 
parameter in 



371 ]. The ratio aj/an {af level density parameter for fission channel, a„ for 
neutron channel) is larger than 1. It decreases towards to an unit with the increase of the 

n 

fissility parameter. The results of a//a„ for a series of nuclei in [20| can be well reproduced. 
Ta is the radius parameter found to be = 1.153 fm [2^. 

With this parametrization 51 fusion-fission reactions have been systematically investi- 
gated with the MWS, Cohen-Swiatecki's, Sierk's and Dahlinger's fission barriers, respec- 
tively, incorporating the proposed approach for describing the fusion (capture) cross sections 
(see Eq.(12)). Calculations of the fission and particle emission widths with the traditional 



statistical theory were introduced in 34!]. The average deviation X\og (see Eq.(16)) of the 
evaporation (and fission) cross sections from the experimental data for these reactions are 
listed in Table II. We find that the average deviation obtained with the MWS potential is 
much smaller than those obtained with the other barriers. By varying the volume, surface 
and curvature coefficients in Eq.(23) and the damping energy E^,, and searching for the min- 
imum of X\og with the MWS fission barriers, we find that the values proposed by Reisdorf 
{20! ] (adopted in the present work) are very close to the corresponding optimal ones. In some 
references the shell damping energy was written as Ed = koA^^^ or similar forms 0, Sol, 40]. 
We find that the minimal deviation is not much improved by changing the value of the 
coefficient /cq- Therefore, in our calculations we consequently keep Reisdorf's coefficients, 
Eq.(23), that contains only one empirically adjustable parameter r^. 
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TABLE II: Average deviation of the evaporation (and fission) cross sections from experimental 
data for 51 fusion-fission reactions with = 1.153 fm. 



model Cohen-Swiatecki Sierk Dahlinger MWS 



? 0.2295 0.2177 0.2373 0.1339 



TABLE III: The minimal average deviation Xio„ and the corresponding optimal value of adopting 



og 

different models for calculating the fission barriers. 



model Cohen-Swiatecki Sierk Dahlinger MWS 
Xiog 0.1813 0.1428 0.1642 0.1086 

r„ 1.106 1.091 1.095 1.120 



To determine the optimal value of r^, we first study the reasonable range of r^. The 
evel density parameter is usually dependent on the nuclear mass number from A/8 to A/12 



g], |38|, |39| . Fig.5 shows the level density parameter a as a function of nuclear mass number A 
adopting different values for ( with Bs = Bk = !)• We estimate the variation region of 
which ranges from about 1.075 to 1.250 fm according to ^4/12 ~ y4/8. Through a variation 
of Ta we can find the optimal values of for a certain model to describe the fission barriers. 
The optimal value of could be different for different fission barrier models. Through 
systematical investigation of the minimal average xfog fusion-fission reactions, 

we search for the optimal parameters set (including the parameters of fission barrier and 
the Ta in level density parameter). The minimal average xfog reactions and the 

corresponding optimal values of for the four fission barrier models are listed in Table III. 
By taking the optimal values of r^, the average deviations xfog froui the experimental data 
get obviously smaller for all of these models, especially for the models of Sierk and Dahlinger. 
The deviation obtained by the modified Woods-Saxon potential is still the smallest one. 

According to the formulas for the fission barrier and the level density parameter, one 
learns that a reasonable calculation of the shell correction S is crucial since the shell correc- 
tion plays a role both for the fission barrier and for the level density parameter, especially for 
heavy systems. Therefore, it is interesting to compare the results with the shell corrections 
obtained by different approaches for searching for the optimal parameters set of the HIVAP 
code. The previous calculations discussed are based on the shell corrections obtained with 
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FIG. 5: (Color online) Level density parameter as a function of nuclear mass number A. 



the 1967 parametrization of Mayers and Swiatecki 3l| for Mld- If we take the shell cor 



rections of MoUer-Nix based on the 1995 parametrization of the macroscopic (liquid-drop) 



energies of nuclei 32|, the minimal deviation Xiog corresponding optimal are 

0.1877 and 1.252 fm with the Sierk's barrier, 0.1681 and 1.268 fm with the MWS fission 
barrier, respectively. Comparing with the results listed in Table III, one finds that using 
the shell corrections based on the 1967 parametrization of the liquid-drop energies Mld of 
nuclei 3l|] the fusion-fission reactions studied in this work can be systematically better de- 
scribed with the present HIVAP code for Wsnv incorporating the proposed approach for ajus- 
Finally, we obtain the optimal parameters set of the HIVAP code: MWS potential model 
for the fission barriers, with = 1.120 fm and together with the 1967 parametrization of 
the liquid-drop energies of nuclei for the shell corrections. 

C. Comparison between the Calculated Results and the Experimental Data 
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With the modified Woods-Saxon potential for the unified description of the entrance 
channel fusion barrier and the macroscopic fission barrier By-^, with Va — 1.120 fm, and 
together with the 1967 parametrization of Mayers and Swiatecki for the shell corrections, 
we obtained the deviations X\og of ^^e evaporation (and fission) cross sections from the 
experimental data for the 51 fusion-fission reactions which are shown in Fig. 6. We find that 
68.3% reactions have values smaller than 0.0714, with which we can estimate the upper and 
lower confidence limits of the systematic errors of the HIVAP code for l^gur (the values are 
1.851^sur and l^sur/1-85, respectively). In the following figures. Fig. 7 - Fig. 12, we present 
the calculated results together with the systematic errors (the shades in the figures) of afus 
and W^sur- The experimental data are also presented for comparison. From these figures, 
one finds that the experimental data can be systematically well reproduced (within about 2 
times deviations) at energies near and above the fusion barriers. 

Fig. 13 shows calculated neutron evaporation residue cross sections for heavy systems 
with ^°^Pb. Because the quasi-fission has not been taken into account in these calculations 
yet, we find that the deviations from the experimental data increase exponentially with 
the increase of Zqn (the positions of the peaks for the evaporation residues can be roughly 
reproduced). This implies that the quasi-fission plays an important role in the reactions 
leading to superheavy nuclei. With the proposed approach for Ucap and M^sur, the ambiguity 
in predicting the probability of quasi-fission could be reduced. 

III. CONCLUSION AND DISCUSSION 

In this work, we proposed a modified Woods-Saxon potential for a unified description of 
the entrance channel fusion barrier and the fission barrier of fusion-fission reactions which 
is based on the Skyrme energy-density functional approach. With the proposed potential 
for the fusion barriers, 120 heavy-ion fusion reactions have been systematically investigated 
together with the barrier penetration concept and an empirical barrier distribution. The 
experimental data for the fusion cross sections afus can be well reproduced and the sys- 
tematic errors are 18% at energies near and above the barriers. Incorporating a statistical 
model HIVAP for describing the decay of the compound nuclei, the evaporation residue 
(and fission) cross sections of 51 fusion- fission reactions have been systematically studied 
simultaneously to investigate and refine some key parameters of the HIVAP code such as 
the fission barrier and the level density parameter. With the optimal value of the radius pa- 
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FIG. 6: (Color online) Deviations xfog °^ calculated evaporation (and fission) cross sections 
from the experimental data for 51 fusion-fission reactions. 

rameter = 1.120 fm of the level density parameter, and with the fission barriers obtained 
by the proposed modified Woods-Saxon potential, the experimental data can be systemati- 
cally reproduced reasonably well. The upper and lower confidence limits of the systematic 
errors of the calculated survival probabilities Wsnr with the HIVAP code are l.SSW^sur and 
VFsur/l-SS, respectively. The infiuence of the shell corrections on the calculated results has 



been explored. The 1967 parametrization of Mayers and Swiatecki [3l| for the macroscopic 
(liquid-drop) energies of nuclei gives better results in the case of these 51 reactions. For the 
systems leading to superheavy nuclei, the infiuence of quasi-fission increases rapidly with 
increasing the compound- nuclear charge number Zcn- With the individual investigation of 
(Jcap and Wsnr, the ambiguity of the prediction of the evaporation cross sections could be 
reduced, which is helpful in testing models for the formation probability Pcn of compound 
nuclei. 

In the present work, the estimated systematic errors based on the enhancement factor 
< < 2 for the capture cross sections at sub-barrier energies are still large, especially 
for heavy systems. A precise prediction of the enhancement factor g and a reduction of the 
corresponding systematic errors are still required, especially for the "cold fusion" in which 
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44| . (TpF denotes the fission cross section, cjer 
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denotes the evaporation residue cross section (a sum over all evaporation channels). The shades in 
this and the following figures denote the systematic errors of the present approach (including both 
the systematic errors of ficap and those of M^sur)) if not otherwise stated. 



the suitable incident energies for producing evaporation residues are near or lower than 
the average fusion barrier. For "hot fusion" systems, the suitable incident energies could 
be higher than the average fusion barrier, and thus the influence of g decreases since the 
capture cross sections are not very sensitive to the enhancement factor g at energies above 
the barrier. In addition, the influence of asymmetric fission and the time- dependent fission 
width 661] have not been taken into account yet. It is known that the nuclear dissipation 
influences the saddle-to-scission time and thus influences the competition between fission 
and particle evaporation. These effects are very important in fission dynamics but they are 
beyond the scope of this work. Work on these aspects is in progress. 
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